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Propeller  Lifting-Surface  Corrections 

By  Wm.  B.  Morgan,1  Member,  Vladimir  Silovic,2  Member,  and  Stephen  B.  Denny,3  Visitor 


‘Correction  factors  for  camber,  ideal  angle  due  to  loading,  and  ideal  angle  due  to 
thickness,  which  ere  based  on  propeller  lifting  surface  theory,  are  presented  for  a 
series  of  propellers.  This  series  consists  of  optimum  free-running  propellers  with 
chordwise  loadings  the  same  as  an  NACA  a  =  0.8  mean  line  and  with  NACA-66 
chordwise  thickness  distributions.  The  results  of  the  calculations  show  that  the  three- 
dimensional  camber  and  ideal  angle  are  generally  greater  than  the  two-dimensional 
camber  and  ideal  angle  at  the  same  lift  coefficient.  The  correction  factors  increase 
with  increasing  expanded  a'  so  ratio,  and  those  for  camber  a  id  ideal  angle  due  to 
loading  decrease  with  increasing  number  of  blades.  Thickness,,  in  general,  induces  a 
positive  angle  to  the  flow,  which  necessitates  a  correction  to  the  blade  pitch.  This 
ideal  angle  is  largest  near  the  blade  root  and  decreases  to  negligible  values  toward 
the  blade  tip  and  increases  with  increasing  number  of  blades.  Skew  induces  an 
inflow  angle,  necessitating  a  pitch  change  which  is  positive  toward  the  blade  root 


and  negative  toward  the  blade  tip. 


Introduction 

The  trend  in  ships  has  been  toward  increasing 
speed,  size,  and  horsepower.  As  a  result,  there 
is  an  increasing  demand  for  the  design  of  propel¬ 
lers  which  are  efficient  and  yet  produce  minimum 
cavitation  and  induced  vibration.  Such  a  task 
requires  a  knowledge  of  the  flow  field  in  which  a 
propeller  operates  and  an  accurate  determination 
of  the  flow  over  the  propeller  blade  surfaces. 

Much  effort  has  been  devoted  recently  to  the 
development  of  a  more  sophisticated  propeller 
theory.  This  theory  has  proceeded  from  the 
simple  momentum  concepts  of  Rankine  [1  ]4 
(1805)  to  the  lifting-line  model  of  Goldstein  (2) 
(1929)  and,  finally,  to  the  lifting-surface  model  of 

1  Naval  Architect,  Hydromechanics  Laboratory,  Naval 
Ship  Research  and  Development  Center,  Washington. 
I).  C. 

*  Naval  Architect,  Hydro-tig  Aerodynamisk  Labora- 
torium,  Lyngby.  Denmark. 

*  Mathematician.  Hydromechanics  Laboratory,  Naval 
Ship  Research  and  Development  Center,  Washington, 
D.  C. 

4  Numbers  in  brackets  designate  References  at  end  of 
paper. 

For  presentation  at  the  Annual  Meeting,  New  York. 
N.  Y.,  November  13-10,  1008,  of  Tiie  Society  of  Navai. 
Architects  and  Marine  Engineers. 


Ludwieg  and  Ginzel  (3)  (1914).  More  recently, 
because  of  the  numerical  evaluations  possible 
with  high-speed  computers,  emphasis  has  been 
placed  on  more  sophisticated  mathematical  models 
than  hitherto  possible;  e.g.,  Sparenberg  [4],  Cox 
[5  J,  Pien  [6],  and  Kerwin  [7j. 

A  complete  review  of  the  various  lifting- 
surface  theories  will  not  be  given  here  since  ade¬ 
quate  reviews  for  work  through  19(54  are  given  by 
Wu  (8),  Isay  [9),  and  Lerbs  et  al.  [10].  Recent 
work  not  mentioned  in  these  three  references  me 
studies  by  Cheng  [II,  12],  Kerwin  and  Leopold 
(13],  Nelson  [14,  15],  Sulmont  (16,  17],  Nishi- 
yama  and  Sasajima  (18],  Malavard  and  Sulmont 
(19),  and  Murray  (20  j.  Murray  develops  the 
lifting-surface  theory  of  both  contrarotating  and 
conventional  propellers. 

In  general,  the  investigators  have  dealt  with  the 
same  basic  mathematical  model  but  have  used 
different  assumptions  to  obtain  numerical  solu¬ 
tions.  The  propeller-blade  boundary  conditions 
are  linearized  and  the  blade  and  its  helical  wake 
are  replaced  by  vortex  systems.  Assumptions 
are  made  for  the  chordwise  and  spanwise  loading 
on  the  blade  and  a  complicated  singular  integral 
is  derived  from  which  the  necessary  distortion 


to  the  blade  can  be  calculated  to  obtain  the  pre¬ 
scribed  loading.  In  addition  to  determining  load¬ 
ing  effects,  Kexwin  and  Leojrald  (21],  Nelson 
[14],  and  Murray  (20)  also  considered  the  effect 
of  blade  thickness,  in  the  linearized  sense,  and 
found  that  it  contributed  significantly  to  the  re¬ 
quired.  blade  distortion. 

For  an  unskewed  propeller  blade,  the  effect  of 
loading  on  the  blade  shape  is  to  require  both  a 
larger  camber  and  a  larger  ideal  angle  than  would 
be  required  in  two-dimensional  flow  in  order  to 
produce  the  same  lift.  The  principal  effect  of 
thickness  is  to  distort  the  flow  such  that  an  in¬ 
crease  in  angle  of  the  blade  is  required  to  maintain 
the  desired  loading.  Likewise,  the  principal 
effect  of  skew  is  to  necessitate  a  blade  angle  change 
but  to  require  little  change  in  camber. 


Because  of  the  linearized  boundary  conditions, 
camber  and  ideal  angle  corrections,  which  are 
independent  of  the  magnitude  of  the  propeller 
loading  [10],  can  be  derived  by  taking  a  ratio  of 
the  three-dimensional  camber  and  ideal  angle 
to  the  two-dimensional  camber  and  ideal  angle, 
respectively.  These  correction  factors  are,  of 
course,  dependent  upon  the  ehordwise  and  span- 
wise  load  (or  pitch)  distributions,  number  of 
blades,  and  blade  area  and  shape.  Similarly,  a 
thickness  correction  factor  can  be  derived  which 
is  independent  of  the  magnitude  of  the  thickness 
but  is,  of  course,  dependent  upon  the  ehordwise 
and  spanwise  thickness  distributions. 

Since  “.ting-surface  correction  factors  can  be 
derived  which  are  independent  of  the  magnitude 
of  the  propeller  loading  and  thickness,  a  system- 


Nomenclafur 


blade  expanded  area 
propeller  disk  area 
blade-thickness  fraction 
coefficient  of  blade  outline 

.  2 rG(r)VD 

lift  coefficient, - — — 

cV, 

two-dimensional  maximum  mean-line  or- 

dmate  for  Cl  =  1.0 

.v.  .  Pd 

power  coefficient, - 

g-l.1V 


thrust  coefficient,  ■ 


-^lolV 


c  *»  section  chord 
P  propeller  diameter 
/  «»  two-dimensional  camber  ordinate 
■=>  maximum  two-dimensional  camber  ordinate 
fp(r,  .Vm„)  «■  maximum  camber  ordinate 
fr(r,  .tv)  *  ehordwise  camber  oidinatc 
fpl(r,  t>)  *  camber  ordinate  induced  by  thickness 
J(r,  tv)  =  ehordwise  camber  normalized  by  maximum 
camber 

G(r)  =>  dimensionless  circulation  from  lifting-line 
theory,  Y(r)/*DV 

G,(r,  8)  =>  dimensionless  circulation  over  lifting  sur¬ 
face.  Ity,  e)/r DV 

G/(r)  *»  dimensionless  circulation  of  helical  free 
vortices 

G /(r,  6}  «•  dimensionless  circulation  of  helical  free 
vortices  over  lifting  surfaces 
k,(r)  =  camber  correction 
k,(r)  =  angle-of-attack  correction  for  thickness 
ka(r)  ideal  angle-of-attack  correction 
»  «=  revolutions  per  unit  time 
Pi>  =  propeller  power 
P,  =  hydrodynamic  pitch 
R  **  propeller  radius 


r  =  radial  coordinate  nondimcnsionalized  by 
propeller  radius 
r*  dimensionless  hub  radius 
(r,8,z)  =  dimensionless  cylindrical  coordinates 
T  =  propeller  thrust 
Hr,  xt)  -  half-thickness  ordinatc- 
lu„(r)  =  maximum-thickness  ordinate 
U(r)  =»  resultant  induced  velocity  from  lifting-line 
theory 

L\(r)  -  velocity  induced  by  loading  normal  to  blade- 
section  chord 

U„j(r)  =  velocity  induced  by  thickness  normal  to 
blade-section  chord 
(Mr)  =  tangential  induced  velocity 
V  =  ship  speed 

I'a(r)  =  speed  of  advance  at  a  given  radius 
V,(P)  =  induced  velocity  vector 
V,(r)  -  resultant  section  inflow  velocity 
w(r)  =  circumferential  mean-wake  coefficient  at  a 
given  radius 

(.v,  y,  s)  =  dimensionless  cartesian  coordinates 

.v,  =  ehordwise  section  abscissa  nondimcnsional¬ 
ized  by  chord 

.vm„  ->  ehordwise  position  of  maximum  camber 
X  =  number  of  blades 
o,(r)  =  section  ideal  angle  of  attack 

=  two-dimensional  ideal  angle  of  attack  for 
Cl  =  1.0 

«,(r)  =  angle-of-attack  correction  from  thickness 
d,(r)  =  hydrodynamic  pitch  angle 

X,  =  induced  advance  coefficient  =  r  tan  d. 

V 

X,  =  apparent  advance  coefficient,  — - 

TtnP 

<r(r,  )  =  source  strength 

ffi(r)  —  angular  position  of  blade-section  leading 
edge 

8.  =  skew  at  blade  tip,  deg 
8t(r)  =  angular  position  of  blade-section  trailing 
edge 
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atic  series  of  these  correction  factors  can  be 
calculated  for  design  use.  By  varying  such 
parameters  as  blade  number,  blade  area,  blade 
skew,  and  hydrodynamic  pitch  angle,  the  effect 
of  these  parameters  for  a  practical  blade  shape 
can  be  ascertained  and  a  better  understanding  of 
the  flow  over  the  blade  will  result.  These  cor¬ 
rection  factors  can  be  utilized  in  propeller  design 
to  obtain  the  blade  pitch  and  camber  without 
empirical  adjustments.  Design  methods — devel¬ 
oped  over  the  years — based  on  empirical  correc¬ 
tion  factors  and  on  correction  factors  from  simpli¬ 
fied  mathematical  models  should  be  regarded  as 
outmoded.5  In  addition  to  providing  the  de¬ 
signer  with  coefficients  with  which  he  can  design 
propellers  by  modern  theoretical  methods,  these 
correction  factors  can  be  used  for  making  qualita¬ 
tive  checks  of  lifting-surface  computer  calcula¬ 
tions,  provided  that  the  prepelier  geometry  and 
loading  do  not  differ  radically  from  the  series. 

Lifting-surface  correction  factors  are  presented 
for  four,  five,  and  six  blades,  for  blade  area  ratios 
from  0.35  to  1.15,  for  hydrodynamic  pitch  ratios 
from  0.4  to  2.0,  and  for  three  blade  skews.  The 
calculations  are  for  optimum"  free-running  pro¬ 
pellers  and  the  results  are  similar  to  those  pre¬ 
sented  by  Lerbs  [10]  but  with  two  essential  dif¬ 
ferences.  One  is  that  the  assumed  chordwise 
load  distribution  is  that  of  the  NACA  «  =  0.-S 
mean  line:  i.e.,  constant  loading  from  the  leading 
edge  to  0.8  chord  and  then  a  constant  slope  to 
zero  at  the  trailing  edge.  The  second  difference 
is  the  inclusion  of  a  correction  for  thickness. 
Experimental  results  have  shown  that  these 
differences  are  significant  in  applying  lifting- 
surface  theory  to  practical  propeller  designs. 
The  thickness  correction  must  be  included  if  the 
propeller  is  to  have  sufficient  pitch,  and  the  NACA 
a  =  0.8  mean  line  is  a  mean  line  which  achieves, 
approximately,  its  theoretical  lift  in  viscous  flow. 
For  example,  the  NACA  a  =  1.0  mean  line 
achieves  only  74  percent  of  its  theoretical  lift 
in  a  real  fluid,  and  it  is  not  feasible  in  the  present 
lifting-surface  theory  to  account  for  viscous 
effects  on  lift.  In  analogy  with  two-dimensional 
results,  the  NACA  a  =  0.8  load  distribution 
requires  an  ideal  angle  for  the  mean  line  to  operate 
at  shock-free  entry. 


1  Reference  (L’L’I  is  an  example  of  an  outmoded  method. 
In  this  particular  reference  the  camber  correction  factors  ki 
and  arc  replaced  by  the  camber  correction  factor  k,. 
and  the  pitch  correction  procedure  is  replaced  by  the  ideal 
angle  due  to  loading  and  thickness. 

•  The  word  optimum  used  here  means  the  propeller 
has  a  constant  hydrodynamic  pitch.  Strictly  speaking, 
such  a  propeller  would  be  optimum  only  when  lightly 
loaded,  when  operating  in  an  inviscid  fluid,  and  if  it  had  an 
infinite  number  of  blades. 


The  following  sections  of  the  paper  will  review 
briefly  the  theory  used  for  the  calculations  and 
outline  the  assumptions  and  limitations.  Results 
of  the  calculations  will  be  discussed.  Some  com¬ 
parisons  between  experiment  and  theory  will 
be  presented  and  conclusions  will  be  drawn  from 
the  apparent  trends  as  to  the  applicability  of  the 
theory. 

Theoretical  Background 
Statement  of  the  Problem 

The  problem  can  be  stated  as  follows:  “For  a 
given  number,  loading,  area,  skew,  chord  distri¬ 
bution,  and  thickness  distribution,  determine 
the  required  blade-section  camber  and  ideal 
angle.”  This  problem  is  similar  to  the  inverse 
problem  of  airfoil  theory,  except  that  the  thick¬ 
ness  is  specified. 

Assumption s 

The  brief  theoretical  background  presented 
here  will  serve  only  to  bring  out  the  salient  fea¬ 
tures  of  the  procedures  used.  Details  of  the 
theories  involved  are  available  in  references 
[12] and  [21]. 

In  the  mathematical  model  for  loading,  a 
distribution  of  bound  vortices  is  assumed  to  cover 
the  blades,  and  free  vortices  are  shed  from  these 
bound  vortices  downstream  along  helical  paths. 
For  thickness,  a  network  of  sources  and  sinks  is 
assumed  to  be  distributed  over  the  blades. 
The  following  assumptions  based  on  this  mathe¬ 
matical  model  are  generally  made: 

1  The  fluid  is  inviscid  and  incompressible. 

2  The  free-stream  velocity  is  axisymmetric 
and  steady,  allowing  the  propeller  to  be  wake- 
adapted. 

3  Each  propeller  blade  is  replaced  by  a  dis¬ 
tribution  of  bound  vortices  for  loading  effects, 
and  sources  and  sinks  for  thickness  effects.  The 
circulation  is  distributed  in  both  the  chordwise 
and  spanwise  directions.  It  follows  from  vortex 
theory  that  free  vortices  are  shed  from  the  bound 
vortices  and,  in  a  coordinate  system  which  ro- 
rates  with  the  propeller,  these  free  vortices  form 
a  general  helical  surface  behind  the  propeller. 

4  Each  of  the  free  vortices  has  a  constant 
diameter  and  a  constant  pitch  in  the  downstream 
direction,  but  the  pitch  may  vary  iu  the  radial 
direction.  This  means  that  effects  of  slipstream 
contraction  and  centrifugal  force  on  the  shape  of 
the  free-vortex  sheets  are  ignored. 

5  The  boundary  conditions  on  the  blade  are 
linearized,  which  implies  that  the  lifting  surface 
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Fig.  1  Propeller  velocity  diagram 


has  only  a  small  deviation  from  the  hydrodynamic 
pitch.  This  assumption  is  similar  to  the  linear¬ 
ized  theory  of  two-dimensional  airfoils  where  the 
boundary  condition  is  not  satisfied  on  the  profile 
but  on  the  profile  chord.  Also,  the  linearization 
enables  separation  of  the  loading  and  thickness 
effects. 

(5  The  pitch  of  the  blade  and  of  the  trailing 
vortex  sheets  is  the  hydrodynamic  pitch  obtained 
from  lifting-line  theory.  Within  the  context 
of  the  linearized  boundary  conditions,  this  is  not 
an  assumption — except  for  skewed  propellers — 
since  the  lifting-line  theory  does  not  account  for 
skew. 

7  The  hub  is  assumed  to  be  small  enough 
that  it  is  not  necessary  to  satisfy  the  hub  boundary 
condition. 

8  Blade  rake  is  not  considered. 

The  assumptions  listed  apply  to  moderately 
loaded  propeller  theory.  Assumption  4  would 
have  to  be  removed  for  heavily  loaded  propellers, 
and  an  additional  assumption  would  have  to  be 
made  for  lightly  loaded  propellers;  namely,  that 
the  effect  of  the  induced  velocities  on  the  pitch  of 
the  vortex  sheets  is  negligible. 

Loading 

Fig.  I  shows  the  velocity  diagram  for  a  pro¬ 
peller  blade  section  in  the  absence  of  thickness. 
From  this  figure  and  within  the  concepts  of  linear¬ 
ized  theory,  the  boundary  condition  at  each  sec¬ 
tion  is 


ai(r)  + 


<>fp(r>  -yc) 
dxc 


- ~  [r,  x<)  ~  yr  (r)  (1) 


where  fp  is  the  camber  along  the  chord  x„  a,  is  the 
ideal  angle  of  attack  (a,  cz  tana,),  V,  is  the  re¬ 


Fig.  2  Coordinate  system 

sultant  inflow  velocity  to  the  blade  section,  U 
is  the  resultant  induced  velocity  from  lifting-line 
theory,  and  UK  is  the  induced  velocity  normal 
to  the  blade-section  chord.  The  angle  /3,  shown 
in  Fig.  1  is  the  hydrodynamic  pitch  angle. 

Both  a  eaK,_r.:an  coordinate  system  (.v,  y,  z)  and 
a  cylindrical  system  (r,  8,  z )  will  be  user!  (see 
Fig.  2).  The  bound  circulation  will  be  assumed 
to  have  the  strength  G,(r)  such  that 

Gr(r,  8)d8  =  G(r)  (2) 

i(r) 

The  coordinates  6,{r)  and  8,(r)  define  the  angu¬ 
lar  position  of  the  blade-section  leading  and  trail¬ 
ing  edges,  respectively,  and  G(r)  is  the  non- 
dimensionalized  bound  circulation  as  determined 
from  lifting-line  theory  (23].  Since  G,(r,  0) 
may  vary  in  both  the  radial  and  chordwise 
directions,  a  free  vortex  may  be  shed  at  each 
point  on  the  blade  surface  and  will  trail  behind 
along  a  helix  of  pitch  jrrtan/3,  (tan  d,  from  lifting¬ 
line  theory). 

It  follows  from  vortex  theory  that  the  strengths 
of  the  helical  free  vortices  are  equal  to  the  negative 
of  the  radial  rate  of  change  of  the  bound  circula¬ 
tion  at  the  point  the  vortex  is  shed.  If  the 
strength  of  each  helical  free  vortex  is  assumed  to 
be  G,{r)  behind  the  trailing  edge,  then 

G,(r)  =  -  ~}r)-  dr  (3) 

dr 

With  this  equation  and  equation  (2),  the  strengths 
of  the  free  vortices  shed  from  the  blade  are  found 
to  be 


4 


Propeller  Lifling-Surfcre  Correction* 


it  f"(,) 

GAr)  =  -  G’;(r,  8)d6dr 

dr 

Or 

X  Y  dr  +  GAr,  6,)  ^  rfr  (4) 
dr  dr 

The  first  term  on  the  right  corrcsj>onds  to  the 
integral  of  the  lsee  vortex  strengths  from  the  lead¬ 
ing  edge  to  the  trailing  edge,  which  results  from  a 
radial  change  in  the  bound  vortices.  The  two 
remaining  terms  on  the  right  are  strengths  of  the 
free  vortices  shed  along  the  blade  outline.  It 
follows  from  this  equation  that,  within  the  lifting 
surface,  the  free  vortices  have  strength 

GAr,  6)  =  -GAr,  0)  ~  dr 
r J,(r>  tv; 

Y  (r,  0o)d6Ar,  0,  >  8  >  0,  (5) 
Je  dr 


From  the  Biot-Savart  law,  the  induced  velocity 
V,  at  any  point  P  on  the  lifting  surface  is  found  to 
be 


W) 

V 


ImiuiiJ  vortiew 


4JL 


free  vortices 

G/(r,0)(S~)  («) 


where  Ai  is  the  area  of  the  lifting  surface,  .4*  is 
the  area  of  the  helical  surface  behind  the  trailing 
edge  of  the  blade,  V  is  the  frecstream  velocity 
or  ship  speed,  S  is  the  vector  distance  fr-u  a 
point  Q  on  the  helical  surface  to  the  point  P 
also  on  the  helical  surface,  d /  is  the  elementary 
vector  tangent  to  the  vortex  line,  and  5  =  |  S  | . 
The  second  term  in  this  equation  can  be  expressed 
in  two  parts  by  equations  (3)  and  (5),  and  the 
quantity 


1  CP  dG(r)  /S  X  d/\ 

2  JJa,  dr  V  S  ) 


dr 


is  now  added  to  and  subtracted  from  this  equation ; 
/la  is  the  area  between  the  trailing  edge  and  a 
generating  line  along  which  a  lifting  line  would  be 
placed.  A  new  coordinate  system  is  introduced 
such  that  the  point  P  always  lies  on  the  lifting 
surface7  ((>].  This  means  that  the  integral 


7  For  this  change  in  coordinate  system,  the  lifting  line 
may  deviate  from  the  blade  surface  at  radial  points  other 
than  /’  when  the  propeller  is  not  of  constant  pitch.  The 
effect  on  the  results  is  probably  small  if  the  pitcli  does  not 
deviate  too  much  from  a  constant  value. 


1  ff 


dG(r)  (SX6I 
dr  \  S3 


M*L 


lifting  line 


is  identical  to  the  velocity  induced  at  a  lifting 
line  by  the  trailing  vortex  system.  The  velocity 
induced  by  the  lifting  surface  itself  is 


(7) 


This  equation  gives  the  velocity  induced  by  a 
single  blade  at  the  point  P.  The  total  velocity 
induced  by  all  the  blades  at  the  point  P  on  one 
of  the  blades  is  obtained  by  summing  this  equa¬ 
tion  with  respect  to  the  blade  position  m,  of  which 
S  is  a  function.  The  velocity  V ,(P)  includes 
the  radial  induced  velocity  as  well  as  the  axial 
and  tangential  velocities.  The  radial  velocity 
does  not  appear  in  the  boundary  condition, 
equation  (1),  and  can  be  neglected.  Then,  the 
normal  velocity  is 

_  V,  =v,  u 

Fulfills  line  V  V, 

and  the  equation  which  must  Ik*  evaluated  in 
determining  bfp/dxc  is 


V,  (  ,  df,\  VAr,  x t)  U(r) 

V  r  +  aj-  >  ■  •  T  ■ 

x(sF)-£T(?rH 


<#> 


where  Z  is  the  number  of  blades.  The  blade 
camber  is  given  by 


fp(r,  x<)  xra,{r) 


.err 

(9) 


To  determine  a,,  this  integral  is  evaluated  from 
the  leading  edge  to  the  trailing  edge  and  fp  is 
taken  to  be  zero. 
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The  equation  for  the  slope  of  the  camber  line, 
equation  (!)),  avoids  integration  limits  of  in¬ 
finity,  but  the  double  integrals  are  difficult  to 
evaluate  because  the  integrands  are  singular 
when  5  becomes  zero.  To  facilitate  the  numerical 
evaluation  of  these  singular  integrals,  Pien  (6] 
and  Cheng  1 12]  made  a  number  of  simplifications: 

1  The  radial  bound  circulation  distribution 
G(r)  and  the  hydrodynamic  pitch  ratio  distribu¬ 
tion  are  expanded  in  a  half-range  cosine  series 
with  a  finite  number  of  terms. 

2  The  singular  points  of  the  integrands  are 
isolated  and  the  integrals  are  evaluated  numeri¬ 
cally.  In  the  region  of  the  singular  point,  the 
radial  integration  is  done  analytically  by  assuming 
that  pitch  and  circulation  are  constant. 

3  For  calculating  the  blade  mean  line,  the 
induced  normal  velocity  is  expressed  as  a  power 
series. 

The  simplifications  are  numerical  approxima¬ 
tions  and  not  assumptions  made  in  the  theory. 
The  precise  accuracy  of  the  evaluations  is  difficult 
to  ascertain. 

Correction  factors  for  camber  and  ideal  angle 
due  to  loading  which  are  independent  of  the 
magnitude  of  propeller  loading  can  be  derived  by 
normalizing  the  three-dimensional  camber  and 
ideal  angle  by  the  two-dimensional  section  values. 
With  regard  to  camber,  only  the  correction  factors 
for  the  maximum  chord  wise  camber  will  be  formed 
in  this  way.  From  equation  (9),  then 

k,(r)  = 

J  max  t'luaxC/j, 


where  .v,„ax  is  the  chordwisc  position  of  maximum 
camber,  CL  is  the  blade-section  lift  coefficient, 
and  /„„t  is  the  two-dimensional  maximum  camber 
and  equal  to 


‘Vmw 

(10) 


fmn\  CmaxC^ 


(11) 


The  function  ClUX  is  the  maximum  mean-line 
ordinate  for  6/  =  1.0  in  two-dimensional  flow. 
For  other  chordwise  positions,  the  three-dimen¬ 
sional  cambers  will  be  normalized  by  the  maxi¬ 
mum  camber;  i.e., 


where  is  the  two-dimensional  ideal  angle  of 
attack  for  C{.  =  1.0  for  the  MAC  A  a  =  0.8  mean 
line. 


Thickness 


Blade  thickness  effects  are  determined  by 
introducing  a  source-sink  system  distributed  over 
the  lifting  surface  [21].  The  induced  velocity 
field  is  given  by 


X  H ,(/>,  r„,  0o)jt~j  dfefo  (14) 
I  wj  1 


where  H,  is  the  velocity  induced  at  the  point  P 
by  a  unit  source  located  at  point  r,  0 ;  i.e. 

H'  *  »?,srad(3s) 

and  where  dxe  is  the  element  of  blade  chord,  and 
<r(r,  6)  is  the  source  strength.  It  is  assumed  that 
the  source  strength  distribution  is  known  and  is 
that  derived  by  the  usual  linear  approximations 
from  airfoil  theory;  i.e. 


<r(ro,  0o)  = 


Vr_  d/(r0,  .vc) 
V  dx(  ~ 


/ra  _  Ut\  1  dl(r0,  x<) 
\X,  K/cosdi  d.ve 


05) 


where  /(r0,  ,\e)  is  the  thickness  distributed  along 
the  chord  ,vt,  Ut  is  the  tangential  induced  velocity, 
and  X,  is  the  apparent  advance  ratio  of  the 
propeller;  i.e.,  V/irnD. 

As  for  the  effect  of  loading,  only  the  axial  and 
tangential  induced  velocity  components  are 
required  from  equation  (14).  If  the  change  in 
camber  along  the  chord  ,vc  due  to  thickness  is 
fpt(r,  xc),  and  the  ideal  angle  induced  by  thickness 
is  a,(r),  the  linearized  boundary  condition  on  the 
blade  is 


«i(r)  + 


£>/p, 


d.vc 


-  (r,  xe) 


(!«) 


/(r,.vc)  = 

JpVi  XtMXj 


(12) 


Correction  factors  for  the  ideal  angle  due  to 
loading  can  be  formed  in  the  same  way  as  the 
camber  correction  factor. 


where  U»,  is  the  induced  velocity  normal  to  the 
section  chord  from  equation  (14)  and  where 
a ,(r)  =  tan  a,(r).  Equation  (16)  is  similar  to 
equation  (1),  and  the  change  in  curvature  due 
to  thickness  is  found  by  integrating  equation 
(1(1)  in  the  manner  of  equation  (9). 
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To  facilitate  evaluating  the  integrals  of  equa¬ 
tion  (M),  Kervvin  and  Leopold  (21  j  used  a  dis¬ 
tribution  of  quadrilateral  source  elements  over  the 
blades.  This  approach  is  similar  to  that  of  Hess 
and  Smith  (24 1  for  potential  flow  about  arbitrary 
bodies,  except  that  the  elements  are  placed  along 
the  chord  rather  than  over  the  surface.  There¬ 
fore,  the  source-sink  distribution  is  not  continuous 
but  is  made  up  of  discrete  elements.  If  a  suffi¬ 
cient  number  of  elements  are  used  to  describe  the 
surface,  the  accuracy  should  be  good. 

Calculations  have  shown  that  thickness  in¬ 
duces  an  inflow  angle  and  a  camber  change  but 
that  camber  alteration  is  small  (21  ( except  for  very 
small  pitch  ratios.  This  flow  distortion  requires 
an  increase  in  blade  angle  (ideal  angle)  to  main¬ 
tain  the  desired  loading.  Correction  factors 
for  this  ideal  angle  can  be  made  independent  of 
the  magnitude  of  thickness  by  dividing  the  angle 
(in  radians)  by  the  blade  thickness  fraction  for 
which  the  calculations  were  made;  e.g. 

<IT> 

where  HTF  is  the  blade  thickness  fraction.  A 
camber  correction  due  to  thickness  can  be  formed 
in  a  similar  manner,  but  since  the  correction  is 
small  it  will  be  ignored.  It  should  be  noted  from 
equation  (17)  that  k,(r )  is  also  a  function  of  the 
propeller  loading  since  V,  is  a  function  of  the 
loading.  This  effect  is  small,  however,  and  no 
correction  will  be  made  in  the  results  which  are 
presented. 

Problems  in  Numerical  Evaluation 

because  of  their  complexity,  the  expressions 
arising  in  the  lifting  surface  calculation  procedure 
described  in  this  paper  were  necessarily  evaluated 
by  digital  computation.  This  resulted  in  re¬ 
peated  numerical  operations  to  obtain  discrete 
values  of  the  functions  involved  and  the  occur¬ 
rence  of  inaccuracies,  which  were  dependent  upon 
the  characteristics  of  the  functions  themselves 
and  were  not  always  predictable  by  the  program 
user. 

The  Cheng  computer  program,  described  in 
reference  (12),  offers  many  input  and  output  data 
options  which  are  helpful  in  describing  propellers 
with  special  geometries  and  loadings.  This 
flexibility,  however,  can  lead  to  inaccuracies  in 
the  calculations.  Points  which  merit  special 
attention  in  both  the  operation  of  the  program 
and  the  analysis  of  the  results  are: 

1  The  computer  program,  as  listed  in  Appen¬ 
dix  If  of  reference  (12),  must  be  modified  slightly 
to  give  correct  results.  The  convergence  criteria 


cited  in  Fortran  statements  201  to  204  of  Sub¬ 
routine  SUB3  should  nad  0.9991)  and  0.0001 
respectively,  rather  than  0.99  and  0.01.  This 
oversight  results  in  considerable  differences  in 
calculated  camber  values  when  compared  to 
results  of  the  former  program  [11]  for  constant 
ehordwise  loading  cases. 

2  Initial  input  to  the  Cheng  program  is  the 
designation  of  a  fan  angle  and  grid  spacing  which 
govern  the  location  and  number  of  control  points. 
The  fan  angle  input  must  contain  the  projected 
view  of  the  propeller  blade,  and  the  grid  spacing 
in  that  fan  should  be  as  fine  as  possible  for  the 
greatest  accuracy  in  flow  calculation  at  the  con¬ 
trol  points.  The  total  number  of  ehordwise 
stations  on  either  side  of  the  control  point,  how¬ 
ever,  cannot  exceed  90.  la  order  to  meet  these 
restrictions,  the  angle  and  grid  spacing  should  be 
carefully  checked  before  each  propeller  calcula¬ 
tion.  and  particular  care  should  be  taken  for 
large  blade  areas  and  highly  skewed  propellers. 
No  automatic  checks  exist  in  the  program  to 
assume  this  task. 

3  Most  critical  of  the  input  parameters  is  the 
specification  of  chord-wise  coordinates  at  which 
the  camber  and  induced  velocities  will  be  cal¬ 
culated.  Considerable  differences  in  calculated 
induced  velocities  appear  for  relatively  small 
abscissa  changes  along  the  chord  near  the  leading 
and  trailing  edges.  The  inclusion  of  these  changes 
in  the  resulting  integrals  yield  rather  small  per¬ 
cent  changes  in  cambei  but  relatively  large  per¬ 
cent  changes  in  ideal  angle.  The  series  dis¬ 
cussed  in  this  paper  was  calculated  with  chord- 
wise  coordinates  specified  from  5.0  to  95.0-percent 
chord.  Questionable  ideal  angle  and  camber 
corrections  arose  occasionally  for  small  blade  areas 
and/or  low  pitch  ratios.  A  change  of  ehordwise 
coordinates  to  a  range  from  2.0  to  98  0  percent 
chord  improved  upon  the  erratic  results  in  these 
regions,  but,  for  consistency  in  the  series  data, 
these  calculated  values  are  not  presented.  It 
should  be  noted  also  that  the  questionable  data 
were  obvious  only  in  comparisons  to  other  data 
in  the  series,  alone  they  perhaps  would  have 
gone  unnoticed.  Independent  propeller  calcula¬ 
tions  should  be  checked  against  the  series  pre¬ 
sented  in  this  paper  or  other  similar  designs  if 
possible. 

4  Lifting-surface  correction  data  arc  given 
herein  only  at  radial  positions  r/R  =  0.3  through 
r/R  =  0.9  for  0.2/f  hub  propellers.  The  lifting- 
surface  calculations  were  made  also  at  r/R  =  0  25 
and  r/R  —  0.95.  However,  the  result0  were 
questionable  which  is  probably  due  to  a  com¬ 
bination  of  numerical  difficulties  involving  the 
singularities  at  the  blade  hub  and  tip. 
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Presentation  e#  the  Series 

Procedure 

The  propellers  for  the  series  were  free-running 
propellers  of  constant  hydrodynamic  pitch  with 
a  hub  diameter  of  0.2  of  the  propeller  diameter. 
Lifting-line  calculations  were  from  a  computer 
program  based  on  the  induction-factor  method 
of  Lerbs  [25].  With  loading  and  pitch  distribu¬ 
tions  available  from  the  lifting-line  calculations, 
lifting-surface  calculations  for  loading  were  made 
using  the  program  developed  by  Cheng  [12]  and 
for  thickness  using  the  program  developed  by 
Kerwin  and  Leopold  [21].  These  have  been 
combined  into  one  program  at  NSRDC  for  run¬ 
ning  on  the  IBM-7090  computer. 

Because  of  the  complicated  nature  of  the  numer¬ 
ical  procedures,  it  is  difficult  to  state  the  accuracy 
of  the  results.  To  this  end,  calculations  were 
made  and  compared  in  Fig.  3  to  those  of  Cox 
[5],  Kerwin  [13],  and  Lerbs  et  a!.  [10]  for  a 
three-bladed  propeller  with  an  induced  advance 
coefficient8  ,\,  of  0.3333,  expanded  area  ratio  of 
0.75,  and  a  constant  chordwise  load  distribution. 
The  camber  correction  factors,  ke(r),  calculated 
using  results  of  the  Cheng  method  are  within  1 
percent  of  those  derived  from  the  Lerbs  method 
for  radii  between  0.4  and  0.8  but  deviate  consider¬ 
ably  at  the  0.3  and  0.9  radii.  Both  the  Cox  and 
the  Kerwin  methods  give  correction  factors  which 
are,  in  general,  higher  than  those  obtained  by  the 
Lerbs  and  the  Cheng  methods.  These  compari¬ 
sons  are  consistent  with  the  way  Cox  and  Kerwin 
arrive  at  their  camber  correction  factors  [10]. 

*  The  induced  advance  coefficient  X,  is  related  to  the 
hydrodynamic  pitch  ratio  by  rX,  =  r  tan$,  =  (P/D),. 


Since  the  calculations  are  based  on  similar  pro¬ 
cedures,  the  reason  for  the  difference  between  the 
results  of  Lerbs  and  Cheng  is  not  known.  The 
difference  toward  the  blade  tip  is  probably  caused 
by  different  numerical  techniques,  and  the  differ¬ 
ence  toward  the  root  is  mainly  due  to  the  lifting¬ 
line  calculations;  i.c.,  induction  factors  versus 
Goldstein  factors. 

For  a  further  comparison  of  lifting-surface 
calculations,  the  camber  correction  factor  for  a 
chordwise  load  distribution  corresponding  to  a 
NACA  a  =  0.8  mean  line  is  also  plotted  in  Fig.  3. 
It  is  interesting  to  nou-  that  the  camber  correction 
factor  for  this  mean  line  is  everywhere  smaller 
than  the  corresponding  factor  for  the  constant 
load  distribution.  The  three-dimensional  camber 
is  higher,  however,  since  Cm**  is  approximately 
23  percent  larger  for  the  a  —  0.8  than  for  the  a  = 
1.0  mean  line. 

Lifting-surface  corrections  were  calculated  for 
propellers  with  four,  five,  and  six  blades,  with 
expanded  area  ratios  from  0.35  to  1.15,  with 
hydrodynamic  pitch  ratios  of  0.4  to  2.0,  and  with 
a  symmetrical  outline  and  skew  angles  of  7 , 
14,  and  21  deg.  The  skew  angle  6,  is  defined  as  the 
angle  between  two  straight  lines  in  the  projected 
plane — one  from  the  shaft  centerline  through 
the  midchord  at  the  root  section  and  the  other 


Table  1  Ordinate*  for  NACA  66  (Mod)  Thickness 
Distribution  and  NACA  a  ~  0.8  Camber  Distribution 


Thickness 

Camber 

Station, 

Ordinate, 

Ordinate 

xr 

t/tm** 

///»** 

0 

0 

0 

0.005 

0.0665 

0  0423 

0.0075 

0.0812 

0.0595 

0.0125 

0.1044 

0,0907 

0.025 

0.1466 

0. 1586 

0.05 

0.2066 

0.2712 

0.075 

0.2525 

0.3657 

0.1 

0.2907 

0.4482 

0.15 

0.3521 

0.5869 

0.2 

0.4000 

0.6993 

0.25 

0.4363 

0.7905 

0.3 

0.4637 

0.8635 

0.35 

0.4832 

0.9202 

0.4 

0.4952 

0.9615 

0.45 

0.5 

0  9881 

0.5 

0.4962 

1  0 

0.55 

0  4846 

0.9971 

0.0 

0.4653 

0.9786 

0.65 

0.4383 

0.9434 

0.7 

0  4035 

0.8892 

0.75 

0.3612 

0.8121 

0.8 

0.3110 

0.7027 

0.85 

0.2532 

0  5125 

0.9 

0.1877 

0.3586 

0.95 

0.1143 

0.1713 

0.975 

0.748 

0.0823 

1.0 

0.0333 

0 
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Table  2  Distribution  of  Blade  Chord 


MHmTosime 


r 

C(r) 

(12 

t  0338 

0.3 

1  8082 

0  4 

1  9048 

0  5 

2.0907 

0.0 

2  1920 

0  7 

2  2320 

0  8 

2  1719 

0.9 

1  8031 

0  95 

1  5302 

1  0 

0 

Table  3  Skew  Distribution  for  14  iNig  Skew  and 
A  of  1.2 


r 

Skew  R 

»  2 

0 

0.3 

0  0037 

0  4 

0  0148 

0.5 

0  0330 

0.0 

0  0004 

0  7 

0.0957 

0.8 

0  1402 

0.9 

0  1949 

1  0 

0.2010 

Table  4  Calculated  Propellers  with  Symmetrical 
Blade  Outline  for  Four,  Five,  and  Six  Blades 


r\i 

0.4 

0.8 

Ar 

*U 

0.35 

0.35 

0.55 

0.55 

0.75 

0.75 

0  95 

0.95 

1 . 15 

1  15 

1.2 

1.6 

2.0 

0.35 

0,35 

0.35 

0.55 

0.55 

0.55 

0.75 

0.75 

0.75 

0.95 

0  95 

0.95 

1.15 

1.15 

!  15 

from  the  shaft  centerline  through  the  tnidchord 
at  the  blade  tip. 


Choice  of  Parameters 

f  Geometric  properties  of  the  propeller  series  were 
chosen  to  be  consistent  with  present-day  practice. 
All  the  propellers  were  designed  to  have  the  same 
chord  wise  load  distribution  as  that  of  the  NACA 
a  =  0.8  mean  line,  i.e.,  a  constant  chordwise  load 
from  the  blade  leading  edge  to  0.8  of  the  chord, 
and  a  constant  slope  to  zero  at  the  trailing  edge. 
This  mean  line  was  chosen  because  viscous  effects 
on  its  lift  are  small  [20],  and  hence  the  potential 
solution  closely  approximates  the  true  physical 
flow.  The  maximum  ordinate  of  this  mean  line 
for  a  Ci.  -  1.0,  C„„x,  is  0.06790,  and  the  two- 
dimensional  ideal  angle  of  attack  for  CL  =  1.0, 
at,„  is  1.5-1  deg  (26).  The  chordwise  thickness 
distribution  chosen  was  the  NACA-G0  section, 
TMB  modification  [27],  which  has  a  desirable 
pressure  distribution.  Table  1  indicates  the 


Fig.  4  Expanded  blade  outlines.  At/ At,  *»  0.75,  Z  «■  5, 
wish  symmetrical  blade  outline  and  14-dcg  skew 


half-thickness  ordinates  as  normalized  by  the 
maximum  thickness.  However,  the  lifting-sur¬ 
face  correction  for  thickness  would  be  expected 
to  be  essentially  independent  of  the  chordwise 
thickness  distribution,  assuming,  of  course,  that 
the  shapes  arc  reasonable. 

The  spanwise  distribution  of  thickness  was 
assumed  to  be  linear  with  respect  to  radius  and 
was  given  by  the  following  equation : 

=  (BTF  -  0.003) (1  -  r)  +  0.003  (18) 

where  BTF  is  the  blade  thickness  fraction.  'I  nis 
choice  was  based  on  an  examination  of  a  number 
of  propellers  of  current  design. 

The  blade  outlines  chosen  were  slightly  wider 
toward  the  tip  than  the  Troost  outline  and  were 
given  by 


r/0-T(j;)  (,9) 

where  C(r)  is  given  in  Table  2  and  (.  1  h:/A  0)  is  the 
expanded  area  ratio.  This  is  a  mathematical 
outline  given  by  Cox  [5]  with  his  constants 
a  =  0.732  and  x  -  0.7,  and  is  essentially  the 
same  outline  as  given  by  Schocnherr  [28]  with  his 
constants  c  =  0.4  and  n  =  0.3.  An  expanded 
blade  for  an  -U/.lo  of  0.75  and  five  blades  is 
shown  in  Fig.  4. 

The  radial  distribution  of  skew  was  chosen  so 
that  the  blade-section  mid-chord  line  follower!  a 
circular  arc  in  the  expanded  plane,  as  shown  in 
Figure  4,  and  was  calculated  from  the  following 
equation: 
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and  is  the  skew  angle.  As  an  example,  skew 
in  the  expanded  plane  as  calculated  from  this 
(20)  equation  is  given  in  Table  3  for  a  skew  angle  of 
14  deg  and  a  n\t  =  1.2. 

Series  of  Propellers  With  Symmetrical  Blade  Outline 

The  propellers  which  were  calculated  for  this 
systematic  series  are  listed  in  Table  4.  Correc¬ 
tion  factors  for  camber,  ideal  angle  due  to  loading, 
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and  ideal  angle  due  to  thickness  as  determined  by 
equations  (10),  (13),  and  (17).  respectively,  are 
listed  in  Table  5  for  these  propellers.  These 
data,  in  general,  are  direct  computer  outputs  and 
arc  not  cross-faired.  The  data  should  not  be 
regarded  as  more  accurate  than  to  the  second 
decimal  place  even  though  three  places  are  given. 
Plots  of  the  correction  factors  versus  the  propeller 


radius  are  shown  in  Figs.  5  through  13  for  a 
representative  part  of  the  data.  Figs.  5,  (i,  and  7 
show  the  correction  factors  for  a  five-bhuled 
propeller  with  a  ttX,  of  1.2  for  the  range  of  ex¬ 
panded  area  ratios  investigated.  Figs.  S,  9,  and 
10  show  the  correction  factors  for  a  five-bladed 
propeller  with  .I*/.lo  =  0.75  for  the  range  of 
induced  advance  coefficients  investigated.  Fi- 
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nally,  Figs.  11,  12,  and  13  show  the  correction 
factors  for  xX,  =  1.2  and  Ak/Aq  for  four,  five, 
and  six  blades. 

In  addition  to  these  results,  the  chordwise 
distributions  of  camber,  as  normalized  by  the 
maximum  camber,  equation  (12),  are  listed  in 
Table  G.  Only  the  normalized  coordinates  for 
part  of  the  series  are  presented  since  the  variation 


with  the  various  parameters  is  small.  The  data 
given  in  this  table  are  direct  computer  output  and, 
obviously,  there  are  discrepancies  at  some  radii. 
These  radii  are  marked  in  the  table.  The 
questionable  data  have  been  carefully  checked 
and  no  apparent  error  could  be  found.  The 
discrepancies  are  probably  linked  to  the  difficulty 
encountered  in  the  calculations  for  low  pitch 
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Table  6  Normalized  Chordwise  Ordinates,  flfr,  xf),  for  Symmetrical  Outline, 
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Table  6  (cont)  Normalized  Chordwise  Ordinates,  Hr,  xj,  for  Symmetrical  Outline, 
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Fig.  5  Camber  correction  factor  for  five  blade;  and 


trailing  edge  of  the  blade,  and  the  differences 
decrease  slightly  toward  the  blade  tip  for  this 
propeller.  As  can  be  seen  from  Table  (5,  there  is  a 
tendency  for  the  coordinates  to  become  flatter 
toward  the  leading  edge  with  increasing  induced 
advance  coefficient  and  with  increasing  expanded 
area  ratio.  Also,  the  ordinates  become  slightly 
fuller  toward  the  trailing  edge  for  increasing  in¬ 
duced  advance  coefficient,  but  show  little  change 
with  area  ratio.  These  data  show  only  slight 
variation  in  the  span  wise  direction,  however, 
and  these  differences  are  essentially  within  the 
accuracy  of  the  calculations.  This  leads  to  the 
interesting  supposition  that  the  normalized  cam- 


Fig.7  Correction  factor  for  ideal  angle  due  to  thickness  ber  ordinates  are  almost  independent  of  their 
for  five  blades,  xX,  -  1.2  and  BTF  =  0.4  spanwise  position. 
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Fig.  10  Correction  factor  for  ideal  angle  due  to  thick¬ 
ness  for  five  blades  and  Ae/Ac  —  0.7  5 


Fig.  11  Camber  correction  factor  for  ir\,  «*  1.2  and 
Ae/Ao  =  0.75 


It  is  quite  apparent  from  Table  5  and  Figs.  5,  (5. 
and  7  that  the  effect  of  expanded  area  ratio 
dominates  the  correction  factors.  The  correc¬ 
tion  factors  all  increase  with  expanded  area  ratio 
for  a  given  number  of  blades  and  pitch  ratio. 
In  general,  the  correction  factors  for  camber  and 


«3  «<  OS  01  0?  Cl  3t 


Fig.  1 2  Correction  ftctor  for  ideal  angle  due  to  loading 
forirX;  »  1.2  and  Ae/Ai  =  0.75 


Fig.  1 }  Correction  factor  for  ideal  angle  due  to  thick¬ 
ness  for  jrX:  =  1.2  and  Al:/Ai,  =  0.75 


ideal  angle  due  to  loading  increase  with  increasing 
induced  advance  coefficient,  except  possibly  near 
the  hub  and  tip  of  the  blades,  but  decrease  for 
increasing  number  of  blades,  i.e.,  for  the  range 
investigated.  The  factors  for  camber  and  ideal 


Table  7  Comparison  of  Chordwise  Distributions  of  Camber  for  Five  Blades, 


7rX|  = 

1 .2,  and  At-/Ao  — 

0.75 

Chord 

Distribution 

Distribution 

Distribution 

Position, 

2- Dimensional 

for 

for 

for 

■Vc 

Distribution 

r  =  0.3 

r  ~  0.6 

r  0.9 

0 

0 

0 

0 

0 

0.025 

0  159 

0  121 

0. 132 

0.133 

0.05 

0  271 

0  230 

0.243 

0  243 

0.1 

0.448 

0  410 

0.424 

0.431 

0.2 

0  699 

0.675 

0  687 

0.698 

0  3 

0  804 

0  856 

0.860 

0  871 

0.4 

0  962 

0,961 

0.963 

0  965 

0  5 

1  (XX) 

1  000 

1  000 

l.(XX) 

0  6 

0  979 

0  977 

0  979 

9  973 

0.7 

0  889 

0  882 

0.889 

0.878 

0  8 

0  703 

0  682 

0.691 

0.682 

0  9 

0.359 

0.364 

0  366 

0.365 

0.93 

0  171 

0  177 

0  177 

0. 176 

0.975 

0  082 

0  085 

0.086 

0  082 

1  0 

0 

0 

0 
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angle  are  greater  than  unity,  except  for  the  higher 
numbers  of  blades  at  low  induced  advance  co¬ 
efficients.  This  means  that  the  three-dimensional 
values  are,  in  general,  greater  than  the  two- 
dimensional  values.  In  most  cases,  the  thickness 
induces  an  ideal  angle  which  is  largest  near  the 
blade  root  and  decreases  to  negligible  values  to¬ 
ward  the  blade  tip.  In  general,  this  angle  de¬ 
creases  with  an  increasing  induced  advance 
coefficient  but  increases  with  increasing  number 
of  blades. 

To  show  the  effect  of  blade  outline,  spanwise 
loading,  and  spanwise  thickness  distribution  on 
the  correction  factors,  propellers  were  calculated 
in  addition  to  those  shown  in  Table  4.  All 
calculations  were  for  a  five-bladed  propeller  with 
an  induced  advance  coefficient  ttX,  of  1.0  and  a 
blade  area  ratio  of  0.75.  Two  blade  outlines,  in 
addition  to  the  one  of  the  series,  were  investigated ; 
one  was  the  Cox  Type  1  and  the  other  was  a 
Troost  outline.  Fig.  14  compares  these  outlines, 
and  Figs.  15-17  show  the  correction  factors  for 
the  two  additional  outlines  and  the  series  outline. 
The  correction  factors  generally  follow  the  blade 
chord  distribution,  i.e.,  the  blade  outline  which  is 
narrowest  near  the  root  and  widest  toward  the 
tip  (Cox  Type  1)  has  the  smallest  correction 
factors  near  the  root  and  the  largest  near  the  tip. 
Correction  factors  for  the  Troost  outline,  which  is 
the  widest  near  the  root  and  the  narrowest  toward 
the  tip,  show  opposite  trends.  This  is  true  for  all 
three  correction  factors,  except  that  the  ideal 
angles  due  to  thickness  are  too  small  near  the 


Fig.  1 5  Comparison  of  camber  correction  factors  for 
three  blade  outlines  of  Fig.  14 


Fig.  16  Comparison  of  correction  factors  for  ideal 
angle  due  to  loading  for  three  blade  outlines  of  Fig.  14 


Fig.  17  Comparison  of  correction  factors  for  ideal 
angle  due  to  thickness  for  three  blade  outlines  of  Fig.  14 


tip  to  show  any  such  effect.  The  trends  of  the 
correction  factors  with  the  blade  outlines  suggest 
that  correction  factors  could  be  approximated 
from  the  series  for  an  arbitrary  outline  by  using 
an  "equivalent”  area  ratio  at  each  radius.  The 
“equivalent”  area  ratio  is  defined  as  the  expanded 
blade  area  ratio  for  a  propeller  of  the  series  which 
has  the  same  chord  at  the  particular  radius  in 
question  as  the  arbitrary  outline.  Correction 
factors  derived  from  the  series  with  "equivalent” 
outlines  for  the  Cox  Type  1  and  Troost  blade 
shapes  were  found  to  be  reasonably  close  to  the 
calculated  values,  but  not  in  every  case  were 
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Fig.  18  Comparison  of  ideal  angle  correction  factors 
for  radial  variation  in  thickness 
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Fig.  19  Comparison  of  three  pitch  distributions  for 
five-bladed  propeller  with  expanded  area  ratio  of  0.75 


they  better  than  the  series  data  without  use  of  an 
“equivalent”  biade  area.  Data  presented  by 
Cox  (5 1  show  similar  trends. 

Calculations  were  also  made  using  a  nonlinear 
thickness  distribution,  and  the  tabulated  thick¬ 
ness  values  for  this  variation  and  the  series  are 
shown  in  Table  X.  The  relation  between  the 
thickness  corrections  for  these  two  distributions  is 
shown  in  Fig.  IS  for  a  five-bladed  propeller  with 
an  induced  advance  coefficient  -X,  of  1.0  and  a 
blade  area  ratio  of  0.75.  A  curve  of  the  correction 
factors  as  estimated  from  the  series  data  is  also 
shown  in  this  figure.  The  estimation  was  made 
by  multiplying  the  series  correction  factor  by  a 
ratio  of  the  local  thickness  given  in  Table  S.  As 
might  be  expected,  the  estimation  gives  a  cor¬ 
rection  which  is  too  low.  However,  if  the  cor¬ 
rection  factors  were  calculated  on  the  basis  of  an 
average  of  the  series  and  the  estimating  pro¬ 
cedure  used,  the  estimated  correction  factors 
would  have  been  extremely  close. 

The  effects  of  two  additional  loading  variations 
were  also  investigated.  These  loading  variations 
resulted  in  a  distribution  of  pitch  which  was  net 
constant.  One  propeller  was  a  wake-adapted 
propeller  and  the  other  had  a  reduced  loading 


Fig.  20  Comparison  of  camber  correction  factors  for 
three  pitch  distributions  of  Fig.  19 


Fig.  21  Comparison  of  correction  factors  for  ideal 
angle  due  to  loading  for  three  pitch  distributions  of 
Fig.  19 


Fig.  22  Comparison  of  correction  factors  for  ideal 
angle  due  to  thickness  for  three  pitch  distributions  of 
Fig.  19 


Table  6  Comparison  of  Radial  Distributions 
of  Thickness 
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(/«,.»  />> 
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0  0215 

o  011139 

n  o 

0  0178 

0  01572 

0  7 

0  0141 

0  01229 

0  8 
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Table  9  Calculated  Propellers  with  Skew 


irX, 

Ae 

- - 0.4- 

0.35  0.75 

1  15 

- - 0.8 - 

0.35  0.75  1.15 

- - 1.2- 

0.35  0.75 

1.15 

- 2.0- 

0.35  0.75 

1.15 

Si.  dr*  | 

7 

7 

7 

7 

7  7 

7 

7  7 

7 

14 

14 

14 

14  14  14 

14 

14 

14 

14 

14 

14 

21 

21 

21 

21 

21 

21 

21 

21 

21 

21 

at  the  tip.  The  resulting  pitch  distributions  are 
shown  in  Fig.  19  along  with  the  wake  distribution 
for  the  wake-adapted  propeller.  The  relation 
between  the  correction  factors  for  these  two  load 
distributions,  as  compared  to  the  series  results, 
is  shown  in  Figs.  20,  21,  and  22.  It  is  quite  ap¬ 
parent  that  the  radial  load  distribution  has  a 
significant  effect  on  the  camber  and  ideal  angle 
due  to  loading  correction  factors.  These  cor¬ 
rection  factors  for  the  propeller  with  the  reduced 
loading  at  the  tip  are  largest  toward  the  blade 
root  and  smallest  toward  the  blade  tip.  The 
correction  factors  for  the  wake-adapted  propeller 
are  close  to  the  series  results,  probably  because 
the  radial  load  distribution  is  much  the  same.  An 
attempt  was  made  to  estimate  the  correction 
factors  from  Figs.  8,  9,  and  10  by  entering  the 
figures  at  the  X(  of  each  radius  for  the  wake- 
adapted  propeller  and  that  with  reduced  loading 
at  the  tip.  Although  the  values  of  the  correction 
factors  estimated  in  this  manner  were,  in  general, 
closer  to  the  calculated  values,  there  was  not  a 
significant  improvement  over  using  the  series  cor¬ 
rection  factors  directly.  The  result  of  the  analysis 
shows  clearly  that  the  effect  of  the  radial  load 
distribution  is  important. 

Series  of  Propellers  With  Skew 

The  propellers  which  were  calculated  with  skew 
are  listed  in  Table  9.  Correction  factors  for 
camber,  ideal  angle  due  to  loading,  and  ideal 
angle  due  to  thickness,  as  determined  by  equations 
(10),  (13),  and  (17),  respectively,  are  listed  in 
Table  10  for  these  propellers.  As  stated  for  the 
symmetrical  blade  outline,  the  data  are  not  cross- 
faired  and  should  not  be  regarded  as  more  accurate 
than  to  the  second  decimal  olace. 

Tables  5  and  10  show  that  the  camber  correc¬ 
tion  factor  k,  is  almost  independent  of  skew  but 
that  it  does  tend  to  increase  slightly  toward  the 
blade  tip  above  the  value  for  the  symmetrical 
outline  for  increasing  skew.  This  effect  increases 
slightly  with  increasing  blade  area  and  with  in¬ 
creasing  induced  advance  coefficient.  The  ideal 
angle  correction  due  to  thickness  tends  to  be 
slightly  smaller  than  the  symmetrical  outline 
value.  Otherwise,  the  relationship  between  this 
correction  factor  and  skew  shows  no  distinct  trend. 


Fig.  23  ideal  angle  correction  factor  induced  by  skew, 
Z  *•  5,irXi  **  1.2,  and  Ag/Ao  **  0.75 


The  dominant  effect  of  skew  is  on  the  ideal  angle 
correction  due  to  loading  ka.  Fig.  23  shows  the 
ideal  angle  correction  factor  induced  by  skew  for 
this  series  versus  the  propeller  radius  for  Z  —  5, 
irXi  =  1.2,  and  A  K//l0  =  0.75.  To  obtain  this 
factor,  the  ideal  angle  correction  factor  for  the 
symmetrical  outline  was  subtracted  from  the 
correction  factor  for  the  skewed  outlines.  This 
figure  shows  that  skew  effect  is  significant  and 
induces  a  positive  angle  toward  the  blade  root 
and  a  negative  angle  toward  the  blade  tip.  This 
effect  has  been  found  by  many  others  (10). 

The  chordwise  distributions  of  camber,  as 
normalized  by  the  maximum  camber,  equation 
(12),  are  listed  in  Table  11  Only  data  for  part  of 
the  series  are  presented,  as  the  results  are  very 
close  to  the  results  obtained  for  the  symmetrical 
outline.  The  trends  shown  by  the  data  for  the 
skewed  outline  are  the  same  as  for  the  symmetrical 
outline. 

Experimental  Cheeks  on  the  Theory 

Using  lifting-surface  theory,  several  propellers 
have  been  designer!  and  tested  at  NSRDC. 
None  of  these  has  been  taken  specifically  from 
the  series,  but  the  computer  programs  used  in 
making  the  lifting-surface  calculations  were  the 
same  as  used  for  the  series  or  were  older  versions 
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Table  1 0  (coni)  Correction  Factors  for  Skewed  Propellers,  Z  =  5 
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Table  1 2  Companion  Between  Theory  and  Experiment  for  a 
Constant  Thrust  Loading 


Propeller 

7. 

•li.AU 

P/D  at 
0.7  R 

DTP 

Cr 

Experi¬ 
mental  X. 
Design  X, 

Experi¬ 
mental  C, 
Design  C , 

A 

5 

1.318 

1.473 

0.054 

0.506 

1.013 

0.970 

B 

3 

0.600 

1.077 

0.010 

a.  563 

0.900 

1.025 

C 

3 

0.600 

1 . 08-1 

0.080 

0.503 

0.987 

0.989 

D 

3 

0.303 

1.080 

0.057 

0.570 

0.995 

1.016 

F 

3 

1.212 

1.073 

0.028 

0.5t0 

1.014 

0.965 

G 

a 

1.480 

1.503 

0.030 

0.301 

1.008 

0.971 

of  the  programs.  Table  12  compares  theoretical 
and  experimental  values  for  six  of  these  propellers, 
designated  A,  B,  C,  D,  F,  and  G;  all  used  the 
NACA  a  =  0.8  mean  line  and  had  no  skew.  In 
general,  the  theoretical  values  are  within  the 
accuracy  of  the  experiments  and  within  the  accu¬ 
racy  with  which  the  blade-section  viscous  drag 
can  be  chosen. 

There  appears  to  be  a  trend  for  wide-bladed 
propellers  A,  F,  and  G  to  !>e  slightly  over-pitched. 
Also,  the  thick  propeller  C  tended  to  be  slightly 
under-pitched  as  compared  to  the  propeller  of 
standard  thickness,  B.  It  is  possible  that  both 
of  these  differences  could  arise  from  determination 
of  the  section  viscous  drag  from  airfoil  data. 
Propellers  A,  F,  and  G  have  blades  with  thickness- 
chord  ratios  lower  than  those  for  which  experi¬ 
mental  data  are  available,  and  propeller  B  has 
blades  with  thickness-chord  ratios  which  are 
higher  than  common.  A  further  discussion  of 
most  of  these  propellers  will  be  found  in  reference 

[29]. 

Conclusions 

From  the  numerous  calculations  made  and  the 
series  data  presented,  a  number  of  conclusions 
can  be  drawn  with  reference  to  lifting-surface 
corrections  for  propellers.  Many  of  these  con¬ 
clusions,  of  course,  have  been  made  previously 
by  other  investigators. 

1  Use  of  a  realistic  chordwise  load  similar  to 
that  for  an  NACA  a  =  0.8  mean  line  results  in 
both  an  induced  camber  and  ideal  angle.  The 
camber  correction  factors  for  this  load  distribution 
are  somewhat  less  than  for  the  constant-load  mean 
line,  but,  in  general,  they  are  greater  than  unity, 
i.e.,  greater  than  the  two-dimensional  camber  for 
the  same  lift.  The  ideal  angle  correction  factor 
due  to  loading  is,  of  course,  zero  for  the  constant 
load. 

2  Correction  factors  can  be  formed  for  the 
camber  and  ideal  angle  due  to  loading  which  arc 
independent  of  the  magnitude  of  propeller  loading. 

3  Thickness  induces  an  angle  and  a  camber, 


but  the  camber  correction  is  negligible  except  for 
low  pitches. 

4  A  correction  factor  can  be  formed  for  the 
ideal  angle  induced  by  thickness  which  is  in¬ 
dependent  of  the  magnitude  of  thickness  but 
dependent  on  the  thickness  distribution. 

5  Correction  factors  for  camber  and  ideal  angle 
due  to  loading  are  largest  near  the  blade  tip  and 
smallest  at  0.4  and  0.5  radii.  They  increase  with 
increasing  expanded  area  ratio  and,  generally, 
with  increasing  induced  advance  coefficient, 
but  they  decrease  with  increasing  number  of 
blades. 

6  The  ideal  angle  due  to  thickness  is  largest 
near  the  blade  root  and  decreases  to  negligible 
values  toward  the  blade  tip.  In  general,  this 
angle  decreases  with  increasing  induced  advance 
coefficient  but  increases  with  increasing  blade 
number. 

7  The  chordwise  distribution  of  camber  is 
somewhat  flatter  toward  the  blade  leading  edge 
and,  in  general,  fuller  toward  the  trailing  edge  as 
compared  with  the  two-dimensional  values.  The 
spanwise  change  in  the  camber  chordwise  dis¬ 
tribution  is  small.  In  practice,  the  use  of  the 
two-dimensional  distribution  is  probably  reason¬ 
able. 

8  The  shape  of  the  blade  outline  has  a  signif¬ 
icant  effect  on  the  correction  factors  for  camber 
and  ideal  angle  due  to  loading. 

9  The  spanwise  load  distribution  has  a  signif¬ 
icant  effect  on  the  correction  factors  for  camber 
and  ideal  angle  due  to  loading. 

10  Skew  has  little  effect  on  the  camber  and 
ideal  angle  due  to  thickness  but  has  a  large  effect 
on  the  ideal  angle  due  to  loading.  Skew  in¬ 
duces  a  positive  angle  near  the  blade  root  and  a 
negative  angle  toward  the  blade  tip. 

1 1  Experimental  results  indicate  that,  for  a 
chordwise  load  distribution  corresponding  to  that 
of  the  NACA  «  =  0.S  mean  line,  the  use  of  lifting- 
surface  corrections  gives  propellers  which,  in 
general,  meet  their  predicted  performance  within 
the  accuracy  of  the  experiments  and  the  accuracy 
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with  which  the  blade-section  viscous  drag  can  be 
chosen. 

12  Correction  factors  derived  herein  should 
replace  those  calculated  by  less  sophisticated 
methods;  for  example,  those  in  reference  (22], 

13  The  numerical  evaluation  of  the  compli¬ 
cated  theoretical  equations  may  lead  to  solutions 
which  are  erroneous,  and  calculations  must  be 
carefully  checked. 
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